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In addition of a flow, plastic deformation of structural glasses (in particular amorphous silica) is 
characterized by a permanent densification. Raman spectroscopic estimators are shown to give a 
full account of the plastic behavior of silica under pressure. While the permanent densification of 
silica has been widely discussed in terms of amorphous-amorphous transition, from a plasticity point 
of view, the evolution of the residual densification with the maximum pressure of a pressure cycle 
can be discussed as a density hardening phenomenon. In the framework of such a mechanical aging 
effect, we propose that the glass structure could be labelled by the maximum pressure experienced 
by the glass and that the saturation of densification could be associated with the densest packing 
of tetrahedra only linked by their vertices. 



I. INTRODUCTION 

The behavior of amorphous silica under high pressure 
has been intensively studied in the last decade*^ 2 -^^^ 
and has recently motivated an increasing amount of nu- 
merical & U (Uea§«2.iML2&imiMMil. Silica appears to 
be elastic up to around 10 GPa and to exhibit a plastic 
behavior at higher pressure. Two features can be empha- 
sized at that level: i) in the elastic regime, the compress- 
ibility exhibits a surprising non monotonous evolution 
with a maximum at around 2-3 GPa^ ' 18 ' 19 ii) when un- 
loading from the plastic regime a permanent densification 
up to 20 % can be observed^ 2 -^^^. No clear evidence 
of alteration of the tetrahedral short range order in this 
unloaded state has been observed**^. 

Above 25 GPa, a change from 4-fold to 6-fold coordi- 
nation is observed^. This 6-fold amorphous seems not to 
be quenchable at zero pressure. When unloading downto 
zero pressure, no trace of 6-fold coordination is obtained. 
Performing X-ray Raman scattering experiments on oxy- 
gen K-edge, Lin et observed a reversible electronic 
bonding transition between 10 and 25 GPa. The lat- 
ter was attributed to a fourfold quartz-like to a sixfold 
stishovitelike change of configuration of silica glass. For 
P > 25 GPa the densification process saturates and af- 
ter unloading to ambient, the density level is the one 
obtained with a maximum pressure P ~ 25 GP a 18 > 19 . 

Questions remain about the nature of the densified 
phase and the mechanism of densification. Recent 



studies^ have proposed the existence of an "activated" 
five-fold coordination at high pressure allowing reorga- 
nization toward a denser tetrahedral network. In for- 
mer studies, in analogy with amorphous ice, Lacksi 2 - 
has proposed a first order transition between 2 differ- 
ent tetrahedral amorphous phases of silica. This tran- 
sition would be kinetically hindered at room tempera- 
ture. This idea of poly-amorphism has received a lot of 
attentionZ«±4J&ik22^121212£. At this stage it is impor- 
tant to separate the known transition at very HP between 
a 4-fold amorphous silica and a 6-fold amorphous silica 
involving a change in the short-range order^ from an 
additional hypothetical transition at lower P between 2 
different amorphous phases of tetrahedral silica and in- 
volving the medium range order—. 

The most recent numerical^ and experimental^ works 
as well as the existence of a continuous range of densities 
for amorphous silica^ after return to ambient pressure 
seem to rule out this idea of transition between 2 different 
forms of amorphous tetrahedral silica. However, the orig- 
inal observation of Lacks remains of interest: perform- 
ing molecular dynamics simulations driven in volume^, 
he noted discontinuities in the pressure signal associated 
with local pressure induced mechanical instabilities. The 
latter are reminiscent of the shear induced mechanical in- 
stabilities previously identified in flowing liquids 2 ^. Note 
that similar localized transitions are widely believed to 
be the main mechanism of shear plasticity of amorphous 
materials 2 ^ and can be associated to the vanishing of 
one eigenvalue of the Hessian matrix of the interatomic 
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potentially. 

In parallel to a structural study, it may be worth con- 
sidering the pressure induced densification process in sil- 
ica according to a mechanical perspective. In particu- 
lar, in the Raman spectroscopic measurements to be pre- 
sented below, we will use cycles of pressure of increasing 
maxima. This protocol will help us to discriminate in 
the spectral patterns, modifications due to the reversible 
elastic deformation of the network from other ones due 
to plastic structural reorganizations. 

We first present the experimental methods i.e. Ra- 
man measurements of silica submitted to pressure cycles 
and analyze their results in terms of mechanical behav- 
ior. Two series of experiments are discussed. In series A 
Raman measurements are performed in situ during the 
loading and unloading stages of successive pressure cy- 
cles. In series B Raman measurements are performed ex 
situ at ambient pressure before and after pressure cycles 
of increasing maximum pressure. 

The results of these spectroscopic measurements are 
presented in the next section and discussed in the frame- 
work of continuum plasticity. The first series of in situ ex- 
periments gives a nice illustration of (densification) plas- 
ticity in the context of a silica glass while the second 
series of experiments allows us to follow the evolution of 
permanent densification vs the pressure maximum of the 
cycle i.e. the density hardening behavior of silica. 

We finally give a discussion, first in terms of mechanical 
behavior, then in terms of amorphous structure. These 
results are of primary importance in the description of 
the mechanical properties of silica. While silica is daily 
used as a calibration sample for nano-indentation mea- 
surements, it appears that the mechanical behavior of 
this material is not fully described. It has been shown 
recently^ that the constitutive laws availabl e 33 ' 34 which 
do not take hardening into account fail to fully describe 
the densification process induced by an indentation test. 
The above data can be used to include hardening in a 
simple constitutive law of silica 3 -^, which gives a precise 
account for this behavior. 

Beyond their interest in terms of mechanical behav- 
ior, these results can also be discussed in terms of amor- 
phous structure. In particular, as detailed below, densi- 
fication can be regarded as a typical glassy phenomenon, 
resulting from a mechanical aging process. Following this 
perspective, the amorphous densified structure would be 
associated with the quench of the structure at high pres- 
sure and could be labelled by a Active pressure in the 
very same acception as a Active temperature can be used 
to label a structure obtained by thermal aging. 



II. EXPERIMENTAL METHODS 

Bars of amorphous silica (Saint-Gobain Quartz IDD) 
are shattered into pieces. Splinters of characteristic 
length of 10 /xm are submitted to cycles of pressure in 
a "Sidoine" Diamond Anvil Cell (DAC) with a maximum 



pressure in the range [1-25 GPa]. Raman spectra are 
collected with a Renishaw RM 1000 micro-spectrometer 
with a Ar + -ion 514 nm, 50 mW laser excitation). A small 
piece of ruby is introduced together with the silica splin- 
ter in order to monitor the pressure level using the shift 
of the Rl-luminescence band. 

Two series of experiments are presented. 

• In series A in situ measurements of the Raman 
spectrum are collected all along the compression 
cycle. Three cycles are presented, The first one 
with P m ax = 7.3 GPa lies in the elastic domain, te 
third cycle consists of a compression up to P max = 
18 GPa followed by a decompression at 1 GPa, fi- 
nally the last cycle consists of a compression up to 
Pmax = 16 GPa and a direct return to ambient 
pressure induced by the breakdown of a diamond. 
Methanol is used as a pressure transmitting fluid in 
this series of experiments (quasi-hydrostatic con- 
ditions). A detailed presentation of these in situ 
experiments can be found in 3 ^. 

• In series B Raman measurements are performed ex 
situ at ambient pressure before and after each cy- 
cle of pressure. A series of compression cycles is 
presented where the pressure maximum is increas- 
ing from 9 GPa for the first cycle to 25 GPa for 
the last cycle. In this case, the pressure transmit- 
ting fluid is a mix similar to 5:1 methanol-ethanol, 
which ensures hydrostaticity up to 16 to 20 GPa 3 -^. 
The Raman spectrum is measured before loading 
and after unloading, the diamond-anvil-cell being 
emptied of the transmitting fluid. 

Figure Q] contrasts the initial spectrum of a sample 
(plain line) and that obtained after a 18 GPa hydrostatic 
loading (red line). 

Between 200 and 750 cm -1 we can identify a main 
band at 440 cm -1 . This band is intense, and affected 
by the densification process: the band gets narrower and 
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FIG. 1: Raman spectra obtained with an amorphous silica 
sample before (plain line) and after (dotted line) a 18 GPa 
hydrostatic loading. 
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is shifted to higher wave numbers. This band has origi- 
nally been attributed to the symmetric stretching mode 
of bridging oxygens between two Si atoms and its evo- 
lution under densification to the decrease of the inter- 
tetrahedral angles Si-O-Si 3 -^. Recent determinations of 
the Raman spectrum from first principles^ have estab- 
lished that this broad band originates more likely from 
bending motions of oxygen atoms that do not belong to 
small rings. 

The defect lines Di and D 2 , at 492 and 605 cm -1 , 
are respectively attributed to the breathing modes of the 
four-membered and three-membered ring a 13 ' 40 ' 41 . Their 
area ratio was previously used in literature as an indica- 
tor for the variation of the "Active temperature " 42 ' 43 ' 44 , 
which is associated with a slight change of density. The 
effect of pressure seems to be better accounted for by the 
shift of the D2 line. As discussed by Polsky et a£^, this 
may be due to pressure induced variations of the Raman 
cross section. The D2 line is of particular interest since 
it has almost no overlap with the main band. Sugiura et 
al.— correlated the position of the D2 line with the ratio 
of the sample density p to its initial density po. The resid- 
ual density evaluated through this relation accounts for 
both irreversible and elastic densification due to residual 
elastic strain. However, several work o 46 ' 47 ' 48 evidenced 
that the D 2 line position is only marginally sensitive to 
residual elastic strains. This correlation was recently 
used to probe the densification gradient surrounding a 
plastic imprint in silica obtained by indentation 3 - 2 -. 



III. RESULTS AND INTERPRETATION 

We describe in the following the experimental results 
and give an interpretation in the framework of the elasto- 
plastic response of continuous media. 

A. Continuum mechanics 

We first recall briefly the formalism of elasto-plasticity. 
Fig. [2] shows the hydrostatic stress as a function of the 
volumetric strain of a medium submitted to pressure. Be- 
low a threshold value of pressure Pq the material remains 
fully elastic and the volumetric deformation is reversible. 
When increasing the pressure beyond threshold, plastic- 
ity sets in and an irreversible deformation adds up to 
the elastic reversible deformation. When unloading from 
Pi > Pq , the material behaves elastically and only the 
elastic part of deformation is recovererd. Loading again 
to a higher pressure P% > P± , we observe the same phe- 
nomenology with a crucial difference: the onset of plas- 
ticity has increased from P Y to . In other words, the 
mechanical behavior depends on the history of mechan- 
ical loading. The material has experienced hardening 4 ^ 
which can be regarded as a mechanical aging. Such a 
behavior is standard for metal shear plasticity and usu- 
ally results at the structural level from the entanglement 
or the pinning of dislocations by impurities^. Metals 
however do not exhibit any volumetric plastic deforma- 




FIG. 2: Sketch of a typical pressure/density curve expected 
in elastoplasticity. A reversible elastic behavior is obtained 
up to the elastic limit pressure Pq . When loading above Pq , 
plasticity sets in and the (elastic) unloading from PY is char- 
acterized by a residual densification. A subsequent loading 
at P2 > P% reproduces the previous unloading curve up to 
P\ before plasticity sets in again. The elastic limit has thus 
evolved under loading from Pq to P 1 and PY . The knowl- 
edge of this density hardening behavior (evolution of the limit 
elastic pressure with density at zero pressure) is necessary to 
give a proper modelling of plasticity of glasses. 
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FIG. 3: Raman shift of the main band ujm in the elastic- 
plastic regime. Successive pressure loading and unloading 
cycles are depicted. Filled and empty symbols correspond 
to loading and unloading curves respectively. For a pressure 
maximum P m ax = 7.3 GPa, a full reversibility is obtained 
(elastic behavior). For a larg er iiictxiiiriiiii pressure Pmax = 18 
GPa, a residual shift is obtained after unloading (plastic be- 
havior). An additional loading at P ma x = 16 GPa reproduces 
the last unloading curve, indicating the increase of the elastic 
limit under loading (hardening). 
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FIG. 4: Raman shift of D2 line ujd2 in the elastic-plastic 
regime (same conditions as above). 



anism). However, irreversible changes of density are fa- 
miliar in soil mechanics and granular materials (dilatancy 
effect)*!. 

The typical mechanical behavior of a material like silica 
glass can be summarized as follows. Before any loading, 
the material is elastic up to a threshold which depends 
both on pressure and shear. Using hydrostatic pressure p 
and equivalent shear stress r as coordinates, this elastic 
threshold corresponds to a continuous curve intersecting 
the two axis. When reaching threshold, if the stress is 
increased, plasticity sets in and the elastic limit is con- 
verted to the maximum value of stress experienced by the 
material. At macroscopic scale, one tries to characterize 
this hardening behavior, relating the evolution of densifi- 
cation with the pressure maximum. At microscopic scale, 
in the present case of silica, in absence of a microscopic 
mechanism as well defined as the motion of dislocation, 
a pending question remains to identify and understand a 
structural signature of density hardening. 



B. Elasto-plastic behavior of silica under pressure 

In this section we present the results of in situ raman 
measurements (series A). Two patterns of the Raman 
spectra measured in situ during the pressure cycles are 
discussed: i) the shift of the main band wm; i) the shift 
of the Z?2 line ud2- 

On Fig. [3] and [4] we show the evolution of the po- 
sition of these two bands for a cycle of pressure up to 
Pmax = 7.3 GPa (black symbols). We observe a full re- 
versibility between loading and unloading. This result is 
consistent with the usual estimate P c ~ 10 GPa for the 
onset of permanent densification in silica. The evolution 
of these indicators is also plotted on the same figures for 
the second series of cycles up to P max = 17.9 GPa. We 
now see that below P c ~ 10 GPa, the same elastic behav- 



FIG. 5: Evolution of the residual densification vs the max- 
imum pressure of the pressure cycle (symbols). The dashed 
line is an indicative sigmoidal curve corresponding to a 18% 
maximum densification. 



ior as before is recovered; then a change of slope on load- 
ing is noticeable at least for the position of the main band; 
then the unloading curve does not reproduce the initial 
loading one and there appears a permanent change of the 
spectrum at ambient pressure; when loading again, one 
follows the very last unloading curve. These two curves 
appear to be very similar to the ideal case of hardening 
plasticity depicted in Fig. [2j This spectroscopic study 
allows us to closely follow the elastic-plastic behavior of 
silica under pressure. 

We now discuss the evolution of permanent densifi- 
cation after cycles of increasing pressure. Even if the 
signal-to-noise ratio is less favorable for the determina- 
tion of the D 2 line than for the main band (see Fig. [3] vs 
Fig. |4]) , we choose the former to estimate the permanent 
densification. As discussed above, the main reason for 
this choice is that in contrast to the main band, the shift 
of the D 2 line appears to be rather independent of the 
elastic stress. This allows us to use it a density probe 
in samples affected by residual elastic stress. Because of 
the microscopic size of the silica samples used in the di- 
amond anvil cell, a quantitative calibration could not be 
performed. Densification can be obtained on millimetric 
samples but require high temperature treatments^ and it 
is far from being obvious that the medium range order 
is directly comparable to the one obtained under high 
pressure at ambient temperature. Following Refi^, the 
densification was estimated using the empirical relation: 



L>2 



^>D2 



(i) 



which was extracted from the experimental data of Sug- 
iura et al.— obtained on shock wave experiments. This 
calibration step is thus only approximative. 

The results are summarized on Fig. [H We thus ob- 
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tain the evolution of the permanent densification with 
the pressure maximum P Y of each pressure cycle. We 
obtain a continuous range of increasing densities with an 
apparent saturation. As shown on Fig. [5J this evolution 
can be approached by a sigmoidal curve. Comparable 
results were obtained recently on window glass using an 
octahedral multi-anvil apparatus^. In mechanical terms, 
this allows us to give a quantitative account of the den- 
sity hardening behavior of silica. In the space of stresses, 
we have obtained the evolution of the maximum pressure 
Py below which the material remains fully elastic as a 
function of the glass density. 

IV. DISCUSSION 

Raman spectroscopic estimators have been shownn to 
give a full account of the density hardening behavior 
of silica under pressure. Raman scattering measure- 
ments during a loading/unloading pressure cycle closely 
reproduce the elastoplastic behavior usually observed in 
stress/strain curves: full reversibility below a limit stress, 
appearance of a residual deformation with an elastic un- 
loading for larger stresses. A crucial difference obviously 
lies in the fact that in the present study the evolution of 
hydrostatic pressure vs density is considered instead of 
shear stress vs shear strain as in usual metal plasticity. 
From a mechanical point of view, the detailed analysis of 
the density increase with the maximum pressure of the 
cycle allowed us to study the evolution of the residual 
densification with the elastic limit pressure P Y . Such a 
knowledge is of crucial importance for the determination 
of constitutive equation modelling the plastic behavior 
of silica. It was shown in Ref.— that when restricting 
the plastic criterion to perfect plasticity (i.e. assuming 
no hardening effect) it was not possible to account for 
the permanent densification of amorphous silica around 
a plastic imprint induced by indentation. Conversely, as 
shown in ReL^ 5 . the assumption of an elliptic plastic crite- 
rion coupling shear stress and pressure together with the 
data of density hardening extracted from the present ex- 
periments allow to describe successfully this phenomenon 
of indentation induced densification. 

From a physical/structural point of view, the present 
results suggest that the description of densified silica need 
not the hypothesis of an amorphous-amorphous tran- 
sition between two types of tetrahedral networks. A 
simpler and alternative scenario consists of pressure in- 
duced reorganizations of the amorphous network allowing 
a more efficient packing of tetrahedra remaining linked 



by their vertices only. Such a scenario does not ex- 
clude the occurence of 5-fold or 6-fold coordinated sil- 
ica in the plastic regime at high pressure. However, the 
latter would correspond to intermediate states between 
two amorphous tetrahedral structures. This occurrence 
of 5-fold or 6-fold coordination would thus simply denote 
the necessity of cutting and rebonding between the two 
structures. The denser structure would thus be quenched 
when pressure decreases down to ambient conditions. In 
that sense, the final structure could be labelled by the 
maximum pressure the material experienced. In this 
context of mechanical aging, the latter pressure could 
be thought of as a "Active" pressure in the same accep- 
tion as the Active temperature in a more classical ther- 
mal aging experiment. Note that these denser structures 
may be affected by internal stresses due to the succes- 
sion of localized reorganizations. More generally, these 
results indicate that the structure and the density of a 
densified sample of vitreous silica will depend crucially 
on the particular path it has followed in the plane pres- 
sure/temperature : it is likely that the medium range 
sructure of densified silica is not fully characterized by 
the only density parameter. 

Looking finally at orders of magnitude for the density 
of the various phases of silica, we observe that the den- 
sity of the stable crystalline equivalent at zero pressure 
(quartz) is 2.65 while the density of the metastable co- 
esite (stable at 2-3 GPa) is 3.01. The density of fused 
silica being 2.2, a 20% increase gives 2.64 i.e values close 
to quartz but well below coesite. In analogy with hexag- 
onal and random close packing, an interesting question 
could thus be whether the maximum observed density of 
amorphous silica corresponds to any geometrical maxi- 
mal packing of tetrahedra bonded by their vertices only 
i.e. the maximum density of the continuous random net- 
work. To our knowledge, though the packing of space 
by tetrahedra, or ellipsoids has been considere d 53 ^ 4 , this 
question has not been discussed yet. 
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